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The statistical properties of the primordial density perturbations has been considered in the past 
decade as a powerful probe of the physical processes taking place in the early universe. Within the 
inflationary paradigm, the properties of the bispectrum are one of the keys that serves to discrimi¬ 
nate among competing scenarios concerning the details of the origin of cosmological perturbations. 

However, all of the scenarios, based on the conventional approach to the so-called “quantum-to- 
classical transition” during inflation, lack the ability to point out the precise physical mechanism 
responsible for generating the inhomogeneity and anisotropy of our universe starting from and ex¬ 
actly homogeneous and isotropic vacuum state associated with the early inflationary regime. In 
past works, we have shown that the proposals involving a spontaneous dynamical reduction of the 
quantum state provide plausible explanations for the birth of said primordial inhomogeneities and 
anisotropies. In the present manuscript we show that, when considering within the context of such 
proposals, the characterization of the spectrum and bispectrum turn out to be quite different from 
those found in the traditional approach, and in particular, some of the statistical features, must be 
treated in a different way leading to some rather different conclusions. 
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I. INTRODUCTION 

Recent advances in observational cosmology are allowing a detailed testing of our understanding of the early universe. 
In particular, the inflationary paradigm, considered as a one of the most important cornerstones of modern cosmology, 
is regarded as providing an explanation for the origin of cosmological structure GH3- Indeed, recent observational 
data (e.g. WMAP [tJ, SDSS f®], Planck 0). is in very good agreement with the theoretical predictions offered by 
the inflationary paradigm. According to this theory, the generation of the primordial inhomogeneities, constituting 
“the seeds of galaxies,” is also rather image-evoking: The perturbations start as quantum fluctuations of the inflaton 
field, and as the universe goes trough an era of accelerated expansion, the physical wavelength associated with the 
perturbations is stretched out reaching cosmological scales. 

At this point one is led to treat the quantum fluctuations as classical density perturbations; there are several kinds 
of arguments for such consideration. However, in most of the original works mentioned above, the issue is only dealt 
implicitly; meanwhile, in more recent works, such as eh mi, specific arguments have been put forward based on the 
highly squeezed nature of the field’s quantum state, or on detailed decoherence arguments jj Finally, the traditional 
picture states that these classical perturbations continue evolving into the cosmic structure responsible for galaxy 
formation, stars, planets and eventually life and human beings. 

However, as has been discussed at length in previous works, the complete theory must not only allow one to find 
expressions that are in agreement with observations, but also, be able to provide an explanation of the precise physical 
mechanism behind its predictions. As originally discussed in [~13j , there is a conceptual difficulty in the standard 
explanation for the birth of cosmic structure provided by inflation: The issue is that according to this account, and 
starting from a highlj0 homogeneous and isotropic state that characterizes the background inflaton and space-time, 
as well as the quantum state of the perturbations, the universe must end in a state with “actual” inhomogeneities and 
anisotropies. Let us be a bit more explicit and recall that, when considering quantum mechanics as a fundamental 
theory applicable, in particular, to the universe as a whole, then one must regard any classical description of the 
state of any system as a sort of imprecise characterization of a complicated quantum mechanical state. The universe 
that we observe today is clearly well described by an inhomogeneous and anisotropic classical state; therefore, such 
description must be considered as an imperfect description of an equally inhomogeneous and anisotropic quantum 
state. Consequently, if we want to consider the inflationary account as providing the physical mechanism for the 
generation of the seeds of structure, such account must contain an explanation for why the quantum state that 
describes our actual universe does not possess the same symmetries as the early quantum state of the universe, which 
happened to be perfectly symmetric (the symmetry being the homogeneity and isotropy). As there is nothing in 
the dynamical evolution (as given by the standard inflationary approach) of the quantum state that can break those 
symmetries, then we are left with an incomplete theory. In fact, this shortcoming has been recognized by others in 
the literature EMI. 

The detailed discussion of the conceptual problems associated with the inflationary paradigm, and the proposed 
paths for resolution of the issue within schemes following standard Quantum Mechanics (e.g., the decoherence program, 
many-worlds interpretation and the consistent histories approach), have been discussed in detail by some of us and 
by others in [H, Oj, [H, [l9| . We will not reproduce those arguments here and invite the interested reader to consult 
those references; the above paragraph is meant only to provide the reader with a brief indication of the kinds of issues 
that a detailed analysis of such questions involves^ 

The idea that has been presented in previous works [12,0, HO-fUl , as a possibility to deal with the aforementioned 
problem, relies on supplementing the standard inflationary model with an hypothesis involving the modification of 
quantum theory so as to include a spontaneous dynamical reduction of the quantum state (sometimes referred as 
the self-induced collapse of the wave function) and to consider such reduction as an actual physical process; taking 
place independently of observers or measuring devices. Regarding the situation at hand, the basic scheme is the 
following: A few e-folds after inflation has started, the universe finds itself in an homogeneous and isotropic quantum 
state, then during the inflationary regime a quantum collapse of the wave function is triggered (by novel physics that 
could possibly be related to quantum gravitational effects), breaking in the process the unitary evolution of quantum 
mechanics and also, in general, the symmetries of the original state. In our situation the post-collapse state will not 
in general be isotropic nor homogeneous. This state is moreover assumed to be sufficiently peaked in the relevant 
variables so as to be susceptible to an approximate classical characterization describing a universe that includes the 
inhomogeneities and anisotropies. That classical description will continue to evolve following the standard physical 
processes of post-inflationary cosmology, into the universe we have today. 


1 We must say that we do not find those arguments to be truly satisfactory, for reasons discussed in detail in fl2| . 

2 At the level of one part in e N where N is the number of e-folds of inflation. 

3 In fact, even if one wants to adopt, say, the many-worlds interpretation, the issue at hand can be rephrased by asking questions about 
the precise nature of the quantum state that can be taken as representing our specific branch of the many worlds. One can also focus 
on the issue of characterizing in a precise mathematical way the quantities that encode the so-called stochastic aspects, which are often 
only vaguely referred to. We will see this in a very concrete way in the the following. 
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The idea of modifying quantum mechanics by incorporating a self-induced collapse is not a new idea [24| and there 
has been a considerably amount of research along this lines: The continuous spontaneous localization (CSL) model 
[ 25 J], representing a continuous version of the Ghirardi-Rimini-Weber (GRW) model [26], and the ideas of Penrose (27} 
and Diosi [28[ regarding gravity as the main agent responsible for the collapse, are among the main programs proposed 
to describe the physical mechanism of a self-induced collapse of the wave function. For more recent examples see Refs. 
3C|. In fact, the implications of applying the CSL model to the inflationary scenario have been studied in Refs. 
32 and 133 leading to interesting results constraining the parameters of the model in terms of the parameters 


29 

31 


characterizing the inflationary model. 

On the other hand, the statistical analyses in Refs. [U, [35( show how the predictions of the simple collapse 
schemes, used in previous works, can be confronted with recent data from the Cosmic Microwave Background (CMB), 
including the 7-yr release of WMAP [36] and the matter power spectrum measured using LRGs by the Sloan Digital 
Sky Survey [37|J. In fact, results from those analyses indicate that while several schemes or “recipes” for the collapse 
are compatible with the observational data, others are not, allowing to establish constraints on the free parameters 
of the schemes. Those works serve to underscore the point that, in addressing the conceptual issues confronting 
the inflationary paradigm, one is not only dealing with “philosophical issues,” but that these have impact in the 
theoretical predictions. While, at the same time, the conclusions drawn can lead to important insights, as well 
as a better understanding of the nature of those predictions and novel ways to consider the relation between such 
predictions and the observations. 

In this work, we will be primarily concerned with the characteristics of the “primordial bispectrum” although we will 
also make some interesting observations regarding the primordial spectrum itself. In the inflationary paradigm, the pri¬ 
mordial bispectrum has been regarded as one of the indicators characterizing possible primordial “non-Gaussianities.” 
By non-Gaussianity, one generically refers to any deviations in the observed fluctuations from the random field of 
linear, Gaussian, curvature perturbations. It is commonly believed that the study of non-Gaussianities will play a 
leading role in furthering our understanding of the physics of the very early universe that created the primordial seeds 
for large-scale structures [H|]. As a matter of fact, the shape and amplitude of the bispectrum is currently being- 
used in order to discriminate among various inflationary models. The amplitude of non-Gaussianities is characterized 
in terms of a dimensionless parameter /nl (defined in Sec. IIVI) . Distinct models of inflation predict different values 
for /nl- The amount of non-Gaussianity from simple inflationary models, which are based on a single scalar filed 
with slow-roll inflationary regime, is expected to be very small j39l - l44| : however, a very large class of more general 
models, such as those with multiple scalar fields, special features in the inflaton potential, non-canonical kinetic terms, 
deviations from Bunch-Davies vacuum, and others (see Ref. [45| and references therein for a review ), are thought to 
be able to generate substantially higher levels of non-GaussianityH 

In this article, we will consider the simplest slow-roll inflation model, within the approach involving the hypothesis 
of “the collapse of the wave function.” We will focus on the simplest collapse scheme involving a single instantaneous 
collapse per mode proposal, although several of the results will be valid in more general schemes. We will also use 
some previous results [4g, [47j, which indicate that certain intrinsic non-linearities associated with the collapse of the 
wave function in a gravitational context, generically induce correlations between the random variables characterizing 
the collapse of the wave function in the different modes. In other words, the collapse of the wave function, cannot 
operate in an exactly mode independent manner, and that whenever the collapse takes place, it leads to correlations in 
the state of the different modes [Z(|. In fact, the basic aspects of the collapse-induced correlation between the modes 
and its particular signature in the CMB angular power spectrum has been analyzed in Ref. [47|]. In the present letter, 
we will be primarily interested in analyze the imprint in the CMB primordial bispectrum, by taking into account the 
correlation between the modes generated by the self-induced collapse. This analysis is made possible because, in our 
approach, we must use explicit variables characterizing the random aspects that determine the primordial bispectrum. 

Therefore, although we agree with the prevailing view, that the primordial bispectrum offers us a valuable obser¬ 
vational tool to examine some physical processes that were relevant in the early universe, the statistical aspects that 
are involved in the comparison of theory and observations, must, in our view, be addressed in a different manner. 

At this point, and in order to avoid any misconceptions, we must warn the reader familiar with the subject that 
we will not be discussing here anything related to the quantum three-point function, but that we are focusing on the 
quantity that the observers actually measure, and which, in our view, differs in essential ways from the former. In 
particular, we will assume a highly Gaussian primordial perturbation, and after including the self-induced collapse 
of the inflaton’s wave function, we will obtain an expression for the observed bispectrum that will be derived from 
the 6th. momentum of distribution associated to the coefficients a; m (corresponding to the harmonic decomposition 
of the temperature anisotropies on the CMB). We note that, in the standard approach, the assumption of an initial 


4 For more details and derivations regarding non-Gaussianity, we refer the reader to the comprehensive review by Komatsu f48tl . Bartolo 
et al. |45ll and others l49H5lH . 
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Gaussian perturbation would lead to a zero value for the bispectrum, up to its cosmic variance. Therefore, for all 
practical purposes, followers of the traditional inflationary paradigm may reinterpret our result as a new estimate 
for the bispectrum variance that is different in shape from the traditional one. In particular, some of the differences 
arises as a result of having included the collapse of the wave function. However from the conceptual point of view, our 
approach and the standard one, differ significantly, hence, what is identified at the theoretical level as the bispectrum 
within our approach is not the same as what is called by that name in the traditional approach. The differences are 
subtle but very important, and we encourage the reader to consult Ref. [l8j |. in order to avoid the misunderstandings 
that result from considering our statistical analysis without appropriately considering the difference in the whole 
approach. 

The reminder of the paper is organized as follows: In Section [TTl we start by reviewing the ideas and technical 
aspects of our proposal, in particular we focus on how to implement the collapse of the wave function during the 
inflationary universe. Then, in Section EU we derive an analytic estimate of the expected value of the observed 
primordial bispectrum within our approach. Afterwards, in Section llVl we discuss the main differences between the 
collapse proposal and the standard approach regarding the estimates of the primordial bispectrum and its statiscal 
aspects. Finally, in Section El we end with a discussion of our conclusions. 

Regarding conventions and notation, we will be using a (—, +, +, +) signature for the space-time metric. The prime 
over the functions f denotes derivatives with respect to the conformal time 77. We will use units where c = h = 1 but 
will keep the gravitational constant G. 


II. BRIEF REVIEW OF THE COLLAPSE PROPOSAL WITHIN THE INFLATIONARY UNIVERSE 

In this section we will present a brief review of the collapse proposal which has been exposed in great detail in 
previous works [H, HD, Sa 53> IIS El] • The main purpose is to present the central ideas behind the proposal in order 
to make the presentation as self- contained as possible. First, we will still describe the formalism in a concise manner; 
although we will not be applying all its details due to its intrinsic complexity, it serves to explain the general idea and 
show explicitly how the collapse generates correlations between different modes. On the other hand, we believe that 
the analysis we will preform, would have analogous counterparts in other much more developed approaches involving 
dynamical collapse theories, such as the GRW [26[ or CSL [25( proposals, and even in approaches that rely on applying 
Bohmian Mechanics to the cosmological problem [s3 . ‘55:]. For examples using CSL in this context see [3ll - l33j . 


A. Semiclassical self-consistent configuration 

Before proceeding with the technical aspects, we wish to discuss the way the gravitational sector and the matter fields 
are treated in our approach. As we have not yet at our disposal a fully workable and satisfactory theory of quantum 
gravity, we will rely on the “semiclassical gravity” approach, which naturally we take only as an effective setting 
rather than something that can be considered as a fundamental theory. The inflationary period is assumed to start at 
energy scales smaller than the Planck mass (~ 10~ 2 — 10 -3 Mp), thus, one can expect that the semiclassical approach 
is a suitable approximation for something that, in principle, ought to be treated in a precise fashion within a quantum 
theory of gravity. The semiclassical framework is characterized by Einstein semiclassical equations G a b = 87rG(T a f,), 
which allow to relate the quantum treatment of the degrees of freedom associated with the matter fields to the 
classical description of gravity in terms of the metric^ The use of such semiclassical picture has two main conceptual 
advantages: 

First, the description and treatment of the metric is always “classical.” As a consequence there is no issue with the 
“quantum-to-classical transition” in the characterization of space-time. Thus we will not need to justify going from 
“metric operators” (e.g. 'L) to classical metric variables (such as ’F). The fact that the space-time remains classical is 
particularly important in the context of models involving dynamical reduction of the wave function, as such “collapse 
or reduction” is regarded as a physical process taking place in time and, therefore, it is clear that a setting allowing 
consideration of full space-time notions is preferred over, say, the “timeless” settings usually encountered in canonical 
approaches to quantum gravity (for some basic references on “the problem of time on quantum gravity” see Ref. ibll ). 

Second, it allows to present a transparent picture of how the inhomogeneities and anisotropies are born from the 
quantum collapse: the initial state of the universe (i.e. the one characterized by a few e-folds after inflation has 
started) is described by the homogeneous and isotropic Bunch-Davies vacuum, and the equally homogeneous and 


5 As wee sill see during the collapse, the semiclassical approximation will not remain 100% valid, because the quantum collapse or jump 
of the quantum state we will have X7 a {T ab ) f 0, while V a G ai> = 0. However, as we will be only interested in the states before and after 
the collapse, this breakdown of the semiclassical approximation would not be important for our present work. 
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isotropic classical Friedmann-Robertson-Walker space-time. Then, at a later stage, the quantum state of the matter 
fields reaches a regime whereby (following the general ideas advocated by R .Penrose) the corresponding state for the 
gravitational degrees of freedom are forbidden, and a quantum collapse of the matter field wave function is triggered 
by some unknown physical mechanism. In this manner, the state resulting from the collapse needs not to share the 
same symmetries as the initial state. After the collapse, the gravitational degrees of freedom are assumed to be, once 
more, accurately described by Einstein semiclassical equation. However, as ( T a b) for the new state needs not to have 
the symmetries of the pre-collapse state, we are led to a geometry that generically will no longer be homogeneous and 
isotropic. 

Another advantage of the semiclassical framework is that the universe can be described by what we have named 
“the semiclassical self-consistent configuration” (SSC), which was originally presented in Ref. (46]. In the following, 
we provide a brief description of such proposal. 

The SSC considers a space-time geometry characterized by a classical space-time metric and a standard quantum 
field theory constructed on that fixed space-time background, together with a particular quantum state in that 
construction, such that the semiclassical Einstein’s equations hold. More precisely, we will say that the set 


{9^{x), <p{x),n(x), JP, |0 e J^} (1) 

corresponds to a SSC if and only if <p(x), tt(x) and Jtf correspond to a quantum field theory constructed over a 
space-time with metric g^(x), and the state |£) in the Hilbert space Jif is such that 


G^[g( x )\ = ^G(^\f liI/ [g(x) > 0(x) > n{x)]\O > ( 2 ) 

for all the points in the space-time manifold. 

The former description is thought to be appropriate in the regime of interests except at those specific times when 
a collapse takes place (because of the reasons mentioned in footnote 4). If we view the SSC as the reasonable 
characterization of situations suitable to be treated within the semi-classical approximation, and consider the various 
post-Planckian cosmological epochs as susceptible to such description, we must consider collapses as taking the 
universe from one such regime say SSC-i to another say SSC-ii. 

The relation between the SSC and the collapse process can be described heuristically in the following way: first, 
within the Hilbert space associated to the given SSC-i, one can consider that a sudden jump |£W) |(A)) target “j s 
about to occur,” with both |£W) and |£W) target j n jfd 1 ). However, in general, the set {, |C^)tar g et} 

will not represent a new SSC. In order to describe a reasonable picture, as presented in Ref. |46| , one needs to relate the 
state | C^)target with another one |C (ii) ) “living” in a new Hilbert space for which {g^, <pW, dT (ii \ |C (ii) )} 

is a valid SSC; the one denoted by SSC-ii. Consequently, one needs to determine first the “target” (non-physical) 
state in to which the initial state is “tempted” to jump, sort of speak, and after that, one can relate such target 
state with a corresponding state in the Hilbert space of the SSC-ii. One then considers that the target state is chosen 
stochastically, guided by the quantum uncertainties of designated field operators, evaluated on the initial state |£W) 
and at the time of collapse. 

In Ref. [46} a prescription for identifying two different SSC’s, during the collapse has been introduced; the idea is 
the following: Assume that the collapse takes place along a Cauchy hypersurface E. A transition from the physical 
state |£W) in to the physical state |C (u ' ) ) in (associated to a certain target non-physical state |C^)tar g et in 

will occur in a way that 


tar g et<C (i) |r«[ff (i) ^ (i) ^ (i) ]IC (i) )tar g etL = 7T (ii) ]|C (ii) > 


( 3 ) 


this is, in such a way that the expectation value of the energy momentum tensor, associated to the states |C^)target 
and |C (u) ) evaluated on the Cauchy hypersurface E, coincides. Note that the left hand side in the expression above 
is meant to be constructed from the elements of the SSC-i (although |C^)tar g et is not really the state of the SSC-i), 
while the right hand side correspond to quantities evaluated using the SSC-ii. 

For the case analyzed in (4fi|, the SSC-i corresponded to the homogeneous and isotropic space-time, i.e. = 0 
and the state of the quantum field corresponding to the Bunch-Davies vacuum, while the SSC-ii corresponded to an 
excitation of a single mode ko, characterized by 41 = F(i )) cos(fco ■ x), with a particular quantum state for the inflaton 
field. The energy momentum tensor associated to that state is compatible with this space-time metric according to 
the SSC recipe. As shown in [46}, the main point is that, when the SSC-ii corresponds to a curvature perturbation 
'F including spatial dependences with wavenumber fco, the normal modes of the corresponding Hilbert space 
which would otherwise be characterized by the typical spatial dependence e lk x , and would be called the k modes, 
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would now contain also corrective contributions of the form e *( fc ±fc 0 ).x s i m p]y because the space-time is not exactly 
homogeneous and isotropic. 

Conversely, this implies that if the energy momentum tensor does not lead to the excitation of other modes besides 
the ±fco in the Newtonian potential 'I' (as required for equation [2] to hold), then, the state of the quantum field needs 
to be excited not just in the modes ±fco, but also in the modes ±nko, with n integer (with the degree of excitation 
decreasing with increasing n). In other words, the post-collapse state is characterized by 

ic 111 )=•• ■ ® id%> ® C> ® icg> ■ ■. w 

This situation can be understood as arising from the intrinsic non-linearity of the SSC construction, which corre¬ 
sponds in a sense to a relativistic version of the Newton-Schrodinger system 1231 ]. 

The details of self-consistent formalism can be consulted in Ref. [46|. We will not use such full fledged formal 
treatment here, in part because the analysis becomes extremely cumbersome even in the treatment of a single mode 
of the inflationary field. Thus, when studying the CMB bispectrum, the task would quickly become a practical 
impossibility. Instead, we will focus on the simpler pragmatical approach first proposed in jl3j but keep in mind the 
lessons learned from the SSC formalism. In the next subsection, we will describe how to implement the basic results 
of the SSC proposal in a practical way. 


B. The primordial curvature perturbation and the self-induced collapse 

We proceed now to introduce the details of the simplest collapse proposal. The starting point will be the same 
as the one underlying the standard slow-roll inflationary model; this is one considers the action of a scalar field (the 
inflaton) minimally coupled to gravity, 

S[<t>, 9ab] = J ~ b<t>g ab - V[0]^ . (5) 

This leads to the Einstein’s field equations G a b = 8nGT a b with T b given by: 


T b a =g ac d c </>cW> + K 


~g cd dc(/)d d (l) - V[cl)] 


( 6 ) 


and the field equations for the scalar field: 


V“V Q 


dvm 

dcj> 


= 0 


(7) 


The next step is to split the metric and the scalar field into a background plus perturbations g a b = g^ b + $9ab, 
(j> = (f> o + 5(j>. The background is represented by a spatially flat FRW space-time with line element ds 2 = a(ii)[—dr] 2 + 
5ijdx' l dx :i ] and the homogeneous part of the scalar field cj>o(r]). 

The scale factor corresponding to the inflationary era is a(r]) ~ —1 /{Hr]) with H the Hubble factor defined as 
H = dt.a/a , thus H ~ const. During inflation H is related to the inflaton potential as H 2 ~ ( 8ttG/3)V. The scalar 
field (f>o (r)) is in the slow-roll regime, which means that <f >g ~ —(a 3 /3a')dcf,V. The slow-roll parameter defined by 
e = ^Mpld^V/V) 2 is considered to be e 1; Mp is the reduced Planck mass defined as Mj, = 1/(8itG). We will 
set a = 1 at the present cosmological time; while we assume that the inflationary period ends at a conformal time 
ry* ~ -HR 22 Mpc. 

Next, we focus on the perturbations. Ignoring the vector and tensorial perturbations, and working in the conformal 
Newtonian gauge, the perturbed space-time is represented by 


ds 2 = a{rf) 2 {— (1 + 2 <&)dg 2 + (1 — 2 ^)Sijdx t dx^], (8) 

with <E> and T functions of the space-time coordinates rj, x l . If we assume that there is no anisotropic stress, Einstein’s 
equations to first order in the perturbations lead to = Vk; thus, combining Einstein’s equations yield: 
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V 2 * = 4tt G^' = -^^6^', (9) 

where in the second equality we used Friedmann’s equations and the definition of the slow-roll parameter. 

Next, we consider the quantization of the theory. As mentioned above, we will work within the collapse-modified 
semiclassical gravity setting. In particular, we will quantize the fluctuation of the inflaton field 5<p(x,r]), but not the 
metric perturbations. For simplicity, we will work with the rescaled field variable y = a5<f>. One then proceeds to 
expand the action ([5]) up to second order in the rescaled variable (i.e. up to second order in the scalar field fluctuations) 


SS (2) = J d 4 x 5 £ ( V 



(Vy) 2 + 



( 10 ) 


The canonical momentum conjugated to y is n = dSC^/dy' = y' — ( a!/a)y = a5(j)'. 

In order to avoid distracting infrared divergences, we set the problem in a finite box of side L. At the end of the 
calculations we can take the continuum limit by taking L —> oo. The field and momentum operators are decomposed 
in plane waves 


y(v, *) = ^E Hr}, Z) = ^ E (ii) 

k k 

where the sum is over the wave vectors k satisfying kiL = 2-Kni for * = 1,2,3 with m integer and y%(rj) = yk{j])d.j: + 
y^{rj)a} - and k^(rj) = gk(v)^k + 9k(v)®' %• The function yk(r]) satisfies the equation: 

y'kiv) + (k 2 - Vk(v) = 0- (12) 

To complete the quantization, we have to specify the mode solutions of m- The canonical commutation relations 
between y and k, will give [ap at ] = L 3 6j: p, when yk{i 7) is chosen to satisfy ykgt ~ Vk9k = * for all k at some time 77. 

The remainder of the choice of 'yk(v) corresponds to the so-called Bunch-Davies (BD) vacuum, which is characterized 
by 


< 13 > 

There is certainly some arbitrariness in selection of a natural vacuum state, but it seems clear that any such natural 
choice would be spatially a homogeneous and isotropic state. The BD vacuum certainly is a homogeneous and isotropic 
state as can be seen by evaluating directly the action of a translation or rotation operator on the state. 

From G a b = 87t G{T a b) and © it follows that 


^ (14) 

It is clear from Eq. (ThO) that if the state of the field is the vacuum state, the metric perturbations vanish, and, thus 
the space-time is homogeneous and isotropic. 

The self-induced collapse model is based on considering that the collapse operates very similar to a kind of self- 
induced “measurement” (evidently, there is no external observer or detector involved). In considering the operators 
used to characterize the post-collapse states, it seems natural therefore to focus on Hermitian operators, which in 
ordinary quantum mechanics are the ones susceptible of direct measurement. We thus separate y^{jf) and kj:{rf) into 
their “real and imaginary parts” yAv) = yj: R (ri) + iyAiv) and %(??) = kj: R (rj) + ik? 1 (ri) . The point is that the 
operators yr’ 1 {v) and k?' 1 (y) are hermitian. Thus@ y = V^-Iyfc^a?' 1 ], E = V / 27£[<jrfc(7?)a?’ 7 ], where 
a R = (ag + a_^)/y/2, al = —*(a£ — a_g)/v / 2- The commutation relations for the a?’ 1 are non-standard 


H[z] denotes the real part of z G C 







( 15 ) 








with all other commutators vanishing. 

Following the a line of thought described above, we assume that the collapse is somehow analogous to an imprecise 
measurement^] of the operators y?' 1 (rf) and ft?’ 1 (jj). The rules according to which the collapse is assumed to happen 
are guided by simplicity and naturalness. 

In particular, as we are taking the view that a collapse effect on a state is analogous to some sort of approximate 
measurement, we will postulate that after the collapse, the expectation values of the field and momentum operators 
in each mode will be related to the uncertainties of the initial state. For the purpose of this work we will work with 
a particular collapse scheme called the Newtonian collapse scheme which is given bj]f| 


(y£’ I m)e = 0, 


(16) 


(7Tj 


R,I 


(v fe))e 





(17) 


where represents the time of collapse for each mode. In the vacuum state, frg is distributed according to a Gaussian 
wave function centered at 0 with spread (Aft^Q. The motivation for choosing such scheme is two-folded. First, the 
calculations performed for this scheme are relatively easier to handle than other schemes that have been considered, 
and second, in Eq. (THl) the variable that is directly related with the Newtonian Potential 'F is the expectation value 
of ft; therefore, it seems natural to consider that the variable affected at the time of collapse is (ftj:(jj |)) while there 
is no change in the mean value of the conjugate variable, thus (y^(r 7 S)) = 0 . 

The random variables x?’ 1 represent values selected randomly from a distribution that in principle might possess 
“small non-Gaussian” features, in the sense that, as discussed in Sec. Ill A1 the collapse induces correlations between 
the modes. In other words, there seems to be a clear justification to presume that the value of the random parameter 
x R ? (with n 7 ^ 1 ) would not be completely independent from the value taken by the random parameters x?’ 1 . 

Heuristically, the point is that once the mode k has collapsed, the modes corresponding to the higher harmonics nk 
become excited as explained above (with the largest effect for n = 2 ), and just as it occurs in the situations involving 
the “stimulated emission of photons” (such as the one underlying the functioning of lasers), a mode that is already 
excited in certain way, has a higher propensity to becoming excited in that manner than it would be otherwise. 
Thus, focusing for simplicity on the modes, which in view of the preceding discussion one expects to be more closely 
correlated, we limit hereafter consideration to the case n = 2 , and characterize the level of correlation by an unknown 
quantity we label A < 1. Specifically, we will assume that the average over possible outcomes of the product of two 
random variables x?’ 1 are characterized by 


x i x E 


- h,%’ + 3 


k,—k f 


A(d 


k,2k' 


k,—2k' 


"2k,k 


’’ + ^-2 


(18a) 


x H = 3 


k.k 1 


-5. 


k, — k' 


A(5 


k,2k' 


-6 


k,—2k' 


J 2k,k> 


-6 


—2 k.k 


>)• 


(18b) 


On the other hand, we must emphasize that our universe corresponds to a single realization of these random 
variables, and thus each of these quantities x?, xl has a single specific value. It is clear that even though we will not 
do that here, one could also investigate how the statistics of x~, xl might be affected in detail by the physical process 


of the collapse. The statistical aspects characterizing these quantities can be studied using as a tool an imaginary 
ensemble of “possible universes,” but as explained in 18] we should, in principle, distinguish those from the statistical 
characterization of such quantities for the particular universe we inhabit; we will discuss these and other aspects in 
more detail the next sections. 


7 An imprecise measurement of an observable is one in which one does not end with an exact eigenstate of that observable, but rather 
with a state that is only peaked around the eigenvalue. Thus, we could consider measuring a certain particle’s position and momentum 
so as to end up with a state that is a wave packet with both position and momentum defined to a limited extent and, which certainly, 
does not entail a conflict with Heisenberg’s uncertainty bound. 

8 In previous works, we have analyzed other collapse schemes such as the independent scheme and the Wigner scheme. See Refs. Bill,HI 
for detailed analyses of the collapse schemes. 
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The next step is to find an expression for the evolution of the expectation values of the field operators at all times. 
This can be done in various ways but the simplest invokes using Ehrenfest’s theorem to obtain the expectation values 
of the field operators at any later time in terms of the expectation values at the time of collapse. The result is: 




cos (fc ?7 - z k ) f l _ l_\ sin(fc ?7 - z k ) / 1 ^ 

k V kr) Zk J k V kijzk 




(19) 


<7rf' r (v))e = ^ cos(fcr; - z k ) + Sm ^ kV Zk ))e, (20) 

with Zk = krjZ. This calculation is explicitly done in Refs. fl3|, [52|. 

Finally, using m, ini) and (1201) we find and expression for the Newtonian potential in terms of the random 
variables and the time of collapse 


( L \ 3/2 ( sm(k<n-z k )\ 

^ (r?) = U ) [ cos ^ ~ ^ +-*- ) 


X 




( 21 ) 


where X g = x5 + ixl, and taking into account Eqs it is clear that 


— ^k.k' + ^k.2k' + 


k,k‘ 


(22a) 


^k^k' ~ X$XZ 2 [<$£ £, + K^k-2k' + ^2fc. —fc' )] ■ 


(22b) 


Expression m is the main result of the present section. It relates the Newtonian potential during inflation to 
the parameters describing the collapse (i.e. the random variables and the time of collapse). It is worth noting 
that all the quantities occurring in (1211) are all complex valued quantities and no quantum operators appear in the 
expression. This is an important difference between what one must confront in our approach and in the standard 
treatment of perturbations during inflation. That is, in the latter approach, the Newtonian potential is strictly a 
quantum operator and then one needs to invoke various kinds of arguments that are often quite vague, and do not 
lead to clearly defined connections with the quantities found in the observations; in particular, one often finds in the 
discussions presented in the standard approaches to the subject, an appeal to quantum randomness that is, however, 
left completely unspecified. Thus, the standard approach suffers from the lack of opportunity for completely clarify 
the characterization regarding the stochastic aspects of the situation (as well as from other conceptual deficiencies that 
have been have discussed in |12|). In our approach, we will not rely on arguments involving horizon-crossing of the 
modes, decoherence or many worlds interpretation of quantum mechanics, to justify the transition from a quantum 
object 'I' to a classical stochastic field ’k, which as we said often leads to rather vague notions about the connection 
of the mathematical expressions used and the objects that emerge from observations. One of the advantages of the 
approach we favor is that, as a result of the collapse postulate, such connection become transparent and specific: it 
is fully encoded in the variables X% characterizing, in a unequivocal manner, all and every stochasticity we will need 
to deal with. 

As is well known, the Newtonian potential is closely related with the the temperature anisotropies whose origins 

can be traced back (in the specific gauge) with the extra red/blue shift photons suffered when emerging from the local 

potential wells/hills. As the values of the two random variables associated to each mode, x5 and xi, are fixed for 

k k 

our universe, it follows from expression eu that these values determine the value of the Newtonian potential Fourier 
components corresponding to our universe, which in turn fix the value of the observed temperature anisotropies. The 
statistic nature of the prescribed distribution of the random variable Xj-_ = x5 + ixi gets transfered to the Newtonian 
potential It is clear that we cannot give a definite prediction for the values that these random variables take in 
our universe, given the intrinsic randomness of the collapse. However, as we will show next, the fact that we have 
a large number of modes k contributing to each of the observed quantities, is the feature that allows a statistical 
analysis and the making of a priori theoretical estimates for the observational quantities. 
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C. Observational quantities 


In this subsection we will present a brief review of the main results obtained in Ref. 1171 which show how the 
correlations between different modes, produced by the self-induced collapse, affect the observational quantities, in 
particular, the CMB angular spectrum. In said work we explored a correlation between the random variables of any 
mode with those of their higher harmonics (which is in a sense reminiscent of the so-called parametric resonances 
found in quantum optics in materials with nonlinear response functions [571 , and just like those has its origins in the 
intrinsic non-linearity of the problem at hand ). As mentioned in j47j |. we found that this effect leads to a departure 
from the standard prediction afflicting mainly the first multipoles of the angular power spectrum. 

The observational quantity of interest corresponds to the temperature fluctuations of the CMB observed today 
on the celestial two-sphere@ The temperature anisotropies are expanded using the spherical harmonics ^(9,<p) = 
Sz, m a i m ^ 5 which means that the coefficients dim can be expressed as 


/ ST 

— {0,ip)Y? m {9, ip)d£l, (23) 

here 9 and <p are the coordinates on the celestial two-sphere, with Yi m (9, ip) the spherical harmonics (l = 0,1, 2 ... and 
—l<m< l), and To ~ 2.725 K the temperature average. 

The different multipole numbers l correspond to different angular scales; low l to large scales and high l to small 
scales. At large angular scales (l < 20), the Sachs-Wolfe effect is essentially the single aspect determining the 
temperature fluctuations in the CMB. That effect relates the anisotropies in the temperature observed today on the 
celestial two-sphere to the inhomogeneities in the last scattering surface, 


^(0,p) ( Vd ,x d ), (24) 

where xd = -Rd( sin 0 sinsinflcos^J, cos 0 ), with Ru the radius of the last scattering surface and tjd is the conformal 
time of decoupling (Rd — 4000 Mpc, rjo — 100 Mpc). The Newtonian potential can be expanded in Fourier modes, 
^ matt (ri D ,x D ) = J2k^ att (VD)e l ^' SD /L 3 . Furthermore, using that e^'* D = ^J2im il jl( kR D)Yi m (6,ip)Y l ^ n (k), ex¬ 
pression (1231) can be rewritten as 


a im = ^^ji(kR D )Yf m (k)*f°- tt (r) D ), (25) 

k 

with ji(kRo ) the spherical Bessel function of order l. 

The Newtonian potential 'F matt appearing in (1251) is evaluated at the time of decoupling which corresponds to the 
matter dominated cosmological epoch. Traditionally, the relation between 4 ,matt and the Newtonian potential at the 
end of inflation is made by making use of the so-called transfer functions T(fc); the transfer functions contain all 
relevant physics from the end of inflation to the latter matter dominated epoch, which includes among others the 
acoustic oscillations of the plasma, and are the source of the main modifications of the simple scale invariant H-Z 
spectra. Thus, 'F5 latt (??£)) = T(k) x Fg, where 4'^ corresponds to the Newtonian potential during inflation and, since 
one is interested in the modes with scales of observational interest, the modes ’F^??) must satisfy —ki] <!C 1 (or, as 
commonly referred, its scale should be “well outside the horizon” during inflation). This is, the coefficients a; m must 
be rewritten as 

A'Tri l ^ 

aim = -^Y'j,( k R D )Y£ n (k)T(k)* s . (26) 

k 

At this point the traditional approach would proceed to calculate averages and higher-correlation functions of the 
coefficients a/ m . Nevertheless, within our model we can make a further step, by substituting Eq. (|21[) (and taking the 


9 The two sphere corresponding to the intersection of our present day past light cone with the last scattering hypersurface 




11 


limit —77 —> 0), which gives an explicit expression for tUg in terms of the parameters of the collapse, in Eq. (1261) . one 
obtainJ^l 


dim — 


£3/2 


'Eg{z k )F hn $)X £ , 


with 


and 


Flm(k) 


ji(kR D )Y l * m (k) T 

fc 3 / 2 


(27) 


(28) 


g{zk) = 


irH F2 

5 M P \ 7 


COS Zk — 


sin Zk 
Zk 


(29) 


Equation m allow us to appreciate one of the advantages of the collapse proposal: the coefficient cq TO , which 
is directly associated with the observational quantities (i.e. the temperature fluctuations), is in turn related to the 
random variables characterizing the collapse. In other words, the statistical features of the coefficients ai m can be 
discussed in terms of the statistics of the random variables Xr. We note that there is no analog expression of 
Eq. (1271) in the standard approach. As a matter of fact, if we follow the conventional way of identifying quantum 
expectation values with classical quantities, the prediction given by the standard inflationary paradigm would be 
(0|^fc|0) = 0 = 'bfc; thus, we would be lead, by Eq. (l26l) . to conclude that ai m = 0; this is, the theoretical prediction 
for the temperature fluctuations would be exactly zero in an evident contradiction 1 j with the observations (see Ref. 
0 ). 


Furthermore, from Eq. (EZD, we can calculate the average over possible outcomes of the random variables of the 
product of two aim coefficients, this is, 


then using Eqs. E^l) 


dirndl'm' 


£3 g(zk)g(zk')Fi m (k)Fi m (k )X^X^,, 
k,k' 


(30) 


9( z k)g(zk')Fi m {k)Fi m {k'){8 % _ % , + A(<%_ 2j p + <%_£,)] 

k,k' 

^2g(z k ) 2 | Fi m (k)F Vrn '{-k ) + \[Fi m (k)Fif m r(-2k) + Fi m (k)Fi> m '(-k/ 2)]| 
k 

= J d 3 k g{z k ) 2 {. F lm (k)F l>ml (k ) + \[F lm {k)F Vrn ,{2k) + F lm {%)F Vm ,(k/ 2)]} , (31) 

where in the last line we used parity of the spherical harmonics, which implies that Fi m {—k) = (—1 ) l Fi m (k) and took 
the limit L —> 00 and k continuous. As we have found in previous works, the data indicates that, to a good degree 
of accuracy (UH^I, the assumption that Zk = kg- is independent of k, i.e. the time of collapse goes as r/S oc /c -1 . 
Therefore, assuming Zk = z = const, and taking T(k) = 1, which as previously indicated is a very good approximation 
for Z < 20, one can compute the integrals in Eq. ©, this is, 


__ 2 i l+l ' 

dlmdl'm' ^”3 

2i l+v 

= 


10 Note that we have multiplied by a factor of 3/(5e) the \I/£ we obtained during inflation, Eq. d21D . This is because, while ^£( 77 ) is 
constant for modes — krj <C 1 during any cosmological epoch, its behavior changes substantially during a change in the equation of state 
for the dominant type of matter in the universe. In particular, during the change from inflation to radiation epochs, '3/ is amplified by 
a factor of approximately 1/e. For a detailed discussion regarding the amplitude within the collapse framework see Ref. [58tl . 

11 Several kinds of arguments would normally be invoked at this point in defense of the standard treatments. For a detailed discussion of 
their merits and shortcomings see Ref. [13. 
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Q'lmQ'l'm' — ( 1) 


with 


Di = 


g (~) 2 1 

(27t) 3 1(1 + 1) 


[1 + AG(/)]; 


(32) 


( 33 ) 


we have also defined G(l) as 


G(l) 


91(1 + 1)VtFT(0 
2 l + 5 / 2 T(l _|_ 3 / 2 ) 


2-^1 (^ — 1 / 2 ; l + 3 / 2 , 1 / 4 ), 


(34) 


where V(l) is the Gamma function and 2 -Fi( a > /3; 7 , z) is the Hypergeometric function F^l 

As we said, a key aspect that in our treatment differs, from those followed in the standard approaches, is the manner 
in which the results from the formalism are connected to observations. This is most clearly exhibited by our result 
regarding the quantity a; m in Eq. m- Despite the fact that we have in principle a close expression for the quantity 
of interest, we cannot use Eq. CZD to make a definite prediction because the expression involves the numbers X? 
that correspond, as we indicated before, to a random choice “made by nature” in the context of the collapse process. 
The way one make predictions is by noting that the sum appearing in Eq. 03 represents a kind of “two-dimensional 
random walk,” i.e. the sum of complex numbers depending on random choices (characterized by the Xg). As is 
well known, for a random walk, one cannot predict the final displacement (which would correspond to the complex 
quantity a; m ), but one might estimate the most likely value of the magnitude of such displacement. Thus, we focus 
precisely on the most likely value of |aj m |, which we denote by |az m |M.L.- In order to compute that quantity, we make 
use of a fiducial (imaginary) ensemble of realizations of the random walk and compute the average over the ensemble 
of the value of the total displacement (\ai m \). Finally we identify: 


|nZm|M.L. — (35) 

Thus, using the previous result, Eq. (l32l) (with l — l’ and m = m 1 ), we find 


|n/ m |M.L. — |nz m p — • (36) 

At this point, one could focus on the quantity that is most often studied in this context, namely 

( 37 ) 

m 

for which we would have the estimate 

(38) 

m 

This is, the angular power spectrum is given by 1(1 + 1)C) M ' L ' = 1(1 + 1 )£);, explicitly: 


1(1 + i)C, M L - = ^ [1 + A G(l)}. (39) 

(27T) a 

It is clear that if A = 0, then the probability distribution functions for the variables X^, Xp are uncorrelated, random, 
and Gaussian (see (I18al) and (I18bl) . Furhtermore, for A = 0 we recover standard flat spectrum (recall that we are 


12 The Hypergeometric function is defined as 2 T’i(“, /3; 7, z) 


r(7) 

r(/3)r( 7 -/3) 


rl t 3-l(l_ t )7-/3-l 

JO (l-tz)° 


dt. 
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Figure 1: Plot of the function G(l) [defined in (l34l) ] in the range 2 < l < 20. The function G(l), which encodes the 
correlation between the modes induced by the collapse, affects only the lowest multipoles 
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Figure 2: Fig. (2a) shows our predicted value for l (l + 1)0 [Eq. (l39l) normalized] in the range 2 < l < 20 considering 
three values of A and the best-fit constant to the observational data from Planck collaboration. In Fig. (2b) we 
show our prediction for 1(1 + 1)0 [Eq. (l39l) normalized] in the range 2 < l < 15 with the best-fit value A = —0.3039 

to Planck data. 


assuming that z = krfi is independent of k), i.e. I(l + 1)0 is independent of l [see Eq. ([39)1 ] . The departure of the flat 
spectrum has a very specific signature that, in principle, can be searched for observationally. The behavior of G(l) 
for 2 < l < 20 is shown in Fig. 1. 

Clearly, by ignoring the transfer function, the traditional prediction of the angular power spectrum is a constant; 
on the other hand, in our approach, the mode correlation induced by the self-induced collapse, altered the standard 
prediction by adding the term A G(l). In Fig. (2a), we show a graphic of our predicted value for 1(1 + 1 )Ci [normalized 
see Eq. (1551) ] with three different values of the A parameter and also a best-fit constant (799.25) to the data from 
Planck collaboration 03 In Fig. (2b) we show a plot for our predicted angular spectrum with the best-fit value 
A = —0.3039 to Planck data. We see that our prediction is, in principle, able to improve the theoretical estimate to 
the observational data. 

As we can see, the effect is stronger for low l and it decreases in a nearly exponential fashion for large l. For a 
more complete analysis, we should take into account the transfer functions T(k ) that are associated with the acoustic 
oscillations originated after the end of the inflationary era, this would lead to the full angular power spectrum. We 
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See http://pla.esac.esa.int/pla/#cosmology 
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want to emphasize that the result in Eq. (1391) is different from the standard approach; first, it contains the function 
g(z), which is a function of the time of collapse, and second, the correlation between different modes provided by the 
collapse is encoded in the function G(l). Evidently, a more detailed analysis should be performed before reaching any 
conclusions, but this modest analysis shows the promising potential of our approach. 

Expression (l39l) is the theoretical estimate to be compared with the observational data, and as should be clear from 
the discussion, the fact that we have to rely on most likely values, for what are in effect the mathematical equivalent 
of random walks, leads us to expect that there should be a general and rough agreement between our estimates and 
observations (assuming the theory is correct). However, we do not really expect a detailed and precise match simply 
due to the intrinsic randomness involved and characterized here in full in terms of the random numbers X. In the 
standard approach, similar considerations involving the randomness of the fluctuations and the uncertainties tied to 
stochasticity, mixed together with arguments involving the limited nature of the region of the universe one is observing, 
also lead to people in the community to expect small differences in predictions and observations. Nevertheless, in 
general, such discussions are based on heuristic arguments; therefore, are limited both in scope and precision. The 
essential difficulty is that, in the standard analysis, the precise stochastic elements are not identified and have no 
mathematical representation in the formalism. We believe that, the fact that in our approach, the stochastic elements 
are clearly identifiable (i.e. the XA, represents a great advantage. For instance, the formalism here provides us with 
an explicit expression for the quantity ai mi such as in Eq. (EH), and, thus, allowing us to study in great detail and in 
a transparent manner, the precise nature of higher order statistical estimates as we will do in the following. 


III. CHARACTERIZING THE PRIMORDIAL CMB BISPECTRUM 

In this section, we will study the connection between the Newtonian potential at the end of inflation and the 
observational quantities obtained from the temperature anisotropies in the CMB; in particular, we will provide the 
connection between the parameters characterizing the collapse and the primordial bispectrum. 

The usual path to look for non-trivial statistical features (e.g. possible non-Gaussianities) in the CMB is to study 
the bispectrum, which is considered to be directly related with the three-point function of the temperature anisotropies 
in harmonic space. The theoretical CMB angular bispectrum is defined as 


- L 'mim2Wi3 


®limi ^l2?ri2 ^137713 • 


(40) 


The bar appearing in (1401) denotes average over an ensemble of universes, but in practice such averaging is not 
viable and is usually “replaced,” in a sense, by averaging over orientations in our own universe; the relation between 
the two types of averages is not clear or direct (this fact has been discussed in great detail in Ref. [HI)- In the 
following, we will show how our approach helps to clarify certain issues that emerge when dealing with the statistical 
aspects of the spectrum and when comparing theoretical estimates and observations. 

Given the definition of the CMB bispectrum and by using considerations involving the assumption of an (in average) 
rotational invariant skvF^I one finds in the literature [Hj, [48| another object called the “angle-averaged bispectrum” 
defined by 


= £ 


h I 2 h 

mi m 2 m 3 


uhhh 

^ 771177127713 


= £ 


h h h 

mi m 2 m 3 


^limi &12 7712 ^13777,3 • 


( 41 ) 


The object 


U I 2 
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m m m J ca ^ ec ^ the Wigner 3 -j symbol (see Ref. [5l| for more details and properties for these 
functions); the bispectrum Bi ± i 2 i 3 is non-vanishing for the values of /,m satisfying the following conditions: 


1. mi + m2 + m3 = 0. 

2 . /1 + h + h is even. 

3. \li — lj | < Ik < h + lj for all permutations of indices. 


14 Here we can see some of the difficulties encountered in standard discussions: If one considers the average over universes, one must face 
not only the fact that this ensemble, even if it is real, is empirically inaccessible, and if it is over orientations then the assumption of 
rotational invariance is simply empty. In fact, if the ensemble is just an imaginary ensemble the rotational invariance can be obtained 
by construction starting with any arbitrary universe and constructing the ensemble by including all possible rotations of the original 
one. It is clear that such considerations have no bearing on the properties of the universe we have access to with our observations. 
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Here it is important to mention that the 2nd. selection rule, comes not from the selection rules corresponding to 
the Wigner 3-j symbols, but from the additional assumption that the 3-point correlation function of the temperature 
anisotropies must be invariant under spatial rotations and translations [591 - 16lj . This is a crucial difference with our 
approach as we will see in the following. Also, the previous three selection rules are called “the triangle conditions” 
as h,l 2 ,h must correspond to the sides of a triangle. As a matter of fact, in the standard approach, one intents to 
estimate Bi 3 i 2 i 3 from the observational data (e.g. see Sec. 3.1 of Ref. ]62j|) by testing different configurations for such 
“triangles.” 

Motivated by the fact that, when following our approach, we can obtain a direct relation between the coefficients 
aim and the random variables characterizing the collapse [Eq. (1271) ] we will be focussing on the expression for the 
“observational” bispectrum: 


B 


obs 

I1I2I3 


E ( h h h 
\ mi m2 m3 

rrii x 


^limi Q'l2‘rri2 ^13 m3 • 


(42) 


Note that this object also contains the Wigner 3-j symbol, therefore, h,l 2 ,h must satisfy the first and the third 
selection rule, meanwhile the 2 nd. rule needs not to be satisfied, since the observed bispectrum is certainly not 
rotational nor translational invariant; the symmetry associated to the 2 nd. selection rule applies to the average not 
to the observed bispectrum which is just one element (realization) of the “ensemble.” As a matter of fact, the more 
general rule, which encompass the 2nd. selection rule, for the Wigner 3-j symbol is h + h + h = integer (or even if 
m 1 = m 2 = m 3 = 0 ). 

The difference between Bi 3 i 2 i 3 and Bf^ is a subtle but important one. While in the definition of B[ 3 i 2 i 3 one 
should perform an average over an ensemble of universes [as is explicitly stated in the definition (1411) ]. the object 
B° 3 l 2 l 3 involves no averages over idealized ensembles whatsoever] 15 ! The only average that is being performed in 
is an average over orientations (i.e. a sum over rrii with a weight given by the Wigner 3-j symbols). 

In our approach, the prediction for B° 3 i 2 i 3 is given explicitly by substituting m in (IT21) which yields 

^ 9 /2 Z. j 9( Z ki)9( z k2)9{ z k3)-^ 1 limi{^l)^l2m2(^‘2)Fl 3rn3 (k3)X^^X^X^. (43) 

ki,k2,k3 


Y 2 obs 


= £ 


h h h 

mi m2 m3 


As is clear from Eq. (143[1 . the collapse bispectrum is in effect a sum of random complex numbers (i.e. a sum where 
each term is characterized by the product X^X^X^, which is itself, a complex random number), leading to what 
can be considered effectively as a two-dimensional (i.e. a complex plane) random walk, in complete analogy with 
what we found in dealing with our estimates of ai m in (1271) . Once more, one cannot give a perfect estimate for the 
final displacement resulting from the random walk. Nevertheless, one might give an estimate of its magnitude. 

Similarly as Bf ^ is characterized by the sum of random variables we cannot give a specific value for its outcome. 
However, as we will see, and in complete analogy of our analysis of the quantities ai m , by focusing on the most likely 
value of the magnitude \B° 3 i s 2 i 3 \ 2 we will obtain a reasonable prediction. 

To recapitulate, the original situation corresponds to the homogeneous and isotropic vacuum state. When a sudden 
change of the initial state takes place due to the collapse (one for each mode), the mode becomes characterized by a 
fixed value of the corresponding random variables; the collection of all the values of such random variables associated 
to all the modes characterizes, therefore, our single and unique universe (so such numbers in consequence completely 
determine \B° 3 i s 2 i 3 \ 2 )', let us denote this set by 


U = {X i ,X i ,. (44) 

Nevertheless, given the stochastic nature of the collapse, we can consider that the universe could have corresponded 
to different set of values for the random variables characterizing the universe in a different manner U = {X^, X^ n ...}. 
The collection of different sets { U , I/,...} thus describe an hypothetical ensemble of universes. We will denote the 
bispectrum associated to the hypothetical ensemble of universes as 
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One should not confuse the fact that when obtaining the specific values of cq m , which result from observations, one needs to perform 
an integral over the CMB sky [as indicated in Eq. (1231 1. with taking averages over ensembles of universes as considered above. 
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BhhhiO = 


E 


h h h 

mi m 2 m 3 


nh+h+h _ 

L9/2 E 5(^i )g(z k 2 )g(z k3 )Fi imi (&)F, 2 m 2 (W 3 m 3 
kl,k 2 ,k 3 


(45) 

with the letter £ labeling the possible set of outcomes corresponding to the random variables that characterizes a 
particular Universe, thus, for example £ = obs, corresponds to the bispectrum observed in our universe, this is Eq. 
da, which we regard as typical member of this hypothetical ensemble. 

Furthermore, we will (just as when dealing with cq m ), make the assumption that the most likely (M.L.) value of 
the magnitude |£>;p 2 z 3 (£)| 2 i R such ensemble is a good estimate of the corresponding one for our own universe, that is 


1 m 


obs 


hh I 


^\B h 


hh 


i(£)Im.L. 


(46) 


Moreover, we can simplify the estimate by taking the ensemble average |-Bi 1 ; 2 i 3 (^)| 2 , i.e. the average over all possible 
outcomes of the bispectrum labeled by £ [Eq. (H5)l ] and identify it with the most likely |Bq hh (01m. L . It is needless 
to say that these two notions are not exactly the same for arbitrary kinds of ensembles; as a matter of fact, the 
relation between the two concepts depends on the probability distribution function (PDF) of the random variables. 
In principle, we do not know the exact PDF, as we have only access to a single realization-our own universe-, but as 
can bee seen from Eqs. m, a very good approximation can be made by assuming that the PDF is highly Gaussian, 
this is, the parameter A which encodes the correlation between the modes caused by the collapse, is small enough that 
a Gaussian PDF for the random variables is justified. Therefore, we can approximate 


I^iW3(0Im. 

which implies that 


(47) 


ra 3 i -i b 'iw 3 (oim. l. * \bim 3 (o\ 2 - 


(48) 


In the reminder of this section, we will focus on computing |'Bi 1 z 2 i 3 (£)| 2 and will omit the label £ with the under¬ 
standing that it represents the ensemble over possible outcomes of £. Therefore, the ensemble average is 


Ty 5 -T 7 \ > / l\ I2 I3 \ I I2 I3 \ -*-*-*- 

I hhh\ 2 -^i l mi m 2 m 3 ) l TO 4 TO 5 7716 ) aiimiai 2 m 2 ai 3 m 3 a hrn 4 a l 2 rn 5 a l 3 rn e - 

m. 1 m.c ' ' ' ' 


(49) 


Let us focus on the quantity a; imi a; 2 m 2 a ; 3 m 3 af im 4 af 2 m 6 af 3m6 , thus using Eq. dUD gives 

a iim 1 a ; 2m2 ai 3 m 3 »* im4 a * 2 m 5 a * 3 m 6 = Je E 9 (z kl )g{z k 2 )g(z k 3 )g(z ki )g(z k 5 )g{z ke ) 

ki,...,k 6 

x Fi imi {k\)Fi 2rn2 (k2)Fi 3 m 3 (k3)F lirni (ki)Fi 2m5 ( &5)F) 3 me (fc6 ) 


X A r Xr Xr Xt Xt Xt . 
k 1 fe 2 k 3 ki kg 


(50) 


As we are assuming that the Xr follow a highly Gaussian PDF, then we can approximate the 6 th momentum 
distribution by: 


Xr Xr Xr Xt X* X* 
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( 51 ) 
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Therefore, substituting Eq. (l5ll) into Eq. (l50l) yields 


a hm 1 ai 2 m2 a i3m 3 a'i irn4 a^ 2rn!i a* 3mB — 

je E 9 (zk 1 )g(z k 2 )g(z k 3 )g(z ki )g(z k 5 )g{z k JFi irni (k 1 )Fi 2 rn 2 (k 2 )F l 3 rn 3 {k 3 )F? imi (k 4 )F l * 2 m; .(k 5 )F l l me (k 6 ) 

ki,...,ke 


x (x^X -£ 2 ■ X^XS ■ X'i Xi + 14 permutations^ 


— ( 'y ' g{z k i)g{z k2 )Fi imi (ki)Fi 2rn2 (k2)X^X^ 3 | | 'y ' g(z k3 )g(zh i )Fi 3m3 [k 3 )F lirn4 (k 4 )X^X. 


J$4 


\Jei ,k 2 


\k 3 ,&4 


x I E 9 {zk 5 )g{z k 6 )F^ mi {k 3 )F^ m 2 {k 6 )X^XZ^ j + 14 permutations 

\k 5 ,k 6 / 

= a hrni ai 2 m 2 ■ ai 3 m 3 a* hmi ■ a([ ams a ^ m6 + 14 permutations, 


(52) 


where the last step uses Eq. ©■ This is an expected result from the beginning because the ai m and the X^ are 
linearly related, hence, the a; m must have the same statistics as the Xg, which we have assumed to be highly Gaussian. 
Thus, all higher order correlators must either be zero (odd) or products of the two-point correlator (even). Of the 15 
permutations, there are 6 of the fo rm aa* ■ aa* ■ aa* and 9 of the form Tlci ■ aa* ■ a*a*. Let us consider the second case 
first and its relation with |H;p 2 j 3 | 2 ; e.g. let us analyze the term 
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^1 ^1 ^3 \ /_ 1 \m 5 

7775 — ?ri5 0 ' 


Al-Dia A 3 ^i,Z 2 ^2,i3(2^1 + f)^3,0? 


(53) 


where in the second line we used the relations given in Eq. (15^1) ; in the fourth line we used the selection rule for 
Wigner 3 -j symbols, and in the last line we used an identity for the sum of the Wigner 3 -j symbols. Consequently, 
since this term is proportional to the monopole (1 = 0), which is unobservable, then it vanishes. The same occurs in 
all 9 terms. Henceforth, all that remains are the 6 terms that we will consider now. Once again, focusing a typical 
term, e.g. 
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= D h D l 2 D l3 , (54) 

where in the second line we used Eq. (1321) and in the last line we used an identity of the Wigner 3 -j symbols. The 
remaining 5 terms are computed in a similar way. The final result is thus 
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\B h i 2l3 \ 2 = D h D l2 D h [1 + S h , l2 (-1) 1 ^ + 5 l2 , h (-l) l '+ 21 * + 5 l3 , h (-l) l *+ 2l ° +26 h , la 5 hlh ] , (55) 

with l = 1,2,.... Consequently, the most likely value for the magnitude of the collapse bispectrum is |/Bq; 2 i 3 |M.L. = 

/-\ 1/2 

[l^hhhl ) > i- e - 


\B h 


I2I3 |M.L. 


J-( H 

r 3 / 2 l lOMpe 1 / 2 


sm z 


cos z - 


[1 + AG , (^i)][l + XG{1 2 )}[1 + AG(l 3 )] 'I ^ 
ll(l\ + l)l2{h + 1 )^ 3^3 + 1 ) / 


x [l + <^Zi,Z 2 (—l) z 3+2ii _|_ d; 2 ; 3 ( —l) il+2i2 + <5j 3i q( —l) ,2+2i3 + 25; 1 i ; 2 i 5/ 2i ; 3 ] 


1/2 


(56) 


where we used the definition of Di and g(z) given in Eqs. (1331) and (0H1) respectively. 

Equation (1551) is the main result of this section. Given the definition of Eq. (H21) . h,l- 2 ,h must correspond to the 
sides of a triangle, otherwise |Bz 1 ; 2 z 3 |m.l. = 0 just as in standard treatments. 

However we must emphasize that \Bi 1 i 2 i 3 \ is not exactly the same as the magnitude of the traditional theoretical 
angle-averaged-bispectrum |Bzp 2 z 3 |, as the latter would correspond to the ensemble average as defined in Eq. (TUI) . In 
fact, in the conventional approach, one would relate an average over an ensemble of universes to a certain quantum 
three-point function, which would of course vanish in the absence of “non-Gaussianities.” In the standard approach 
argument then proceeds, by identifying the average over the ensemble of possible universes with the angle-averaged- 
bispectrum Bi x i 2 i 3 , which is just a orientation average in the only universe we have access to. Thus the theoretical 
quantity is connected with suitable averages over m’s of the quantity aj imi o; 2 TTl 2 a ; 2m3 measured in our own (single) 
universe. 

Our point is, thus, connected with the fact that_ these series of identifications do not have a fully transparent and 
clear justification (as discussed in detail in Ref. [H|); and this, in our view, makes it very difficult to deal, in a clear 
way, with existing and potential delicate issues. One example of such difficulty arises when wanting to study issues 
concerning say large deviations from isotropy in the observed CMB. On the other hand, within the collapse model, we 
find a prediction for the most likely value of |/Bzp 2 z 3 |m.l. which can be related directly with the actual observational 
quantity: 


I B 


llhh 


Actual obs — 


E f h h h 

\ mi 777-2 7773 

rrii ' 



irai ®Z2?ri2 ^ 2^-3 / Actual obs j 


(57) 


in a much more direct and transparent manner. Note, however, that we make a distinction between the theory’s 
prediction for the most likely value of the observational quantity I-Bzp^Jm.l. and the actually observed quantity itself 
|-BziZ 2 z 3 (Actual obs- This is exactly analogous to the distinction beteen a theory prediction of the most likely value of 
the total displacement in a random walk, and the actual displacement observed in a single realization of said walk. 

At this point it is also convenient to clarify certain aspects within our approach. Some readers may claim that our 
expression for the observed bispectrum is equivalent to the cosmic variance of the bispectrum, within the standard 
approach, in the “weak non-Gaussian” limit. In other words, the fact that we used the approximation in Eq. m, 
corresponding to a highly Gaussian PDF for the random parameters of the collapse, could lead some people to state 
that what we are computing is actually the variance of the traditional bispectrum, Eq. SID. when = 0. This 

is partially correct and a few remarks are in order. 

If A = 0, which means that the PDF of the random collapse parameters is exactly Gaussian, and if h + I 2 + h 
is even, then our expression for \ Buu u Im t,. coincides exactly with the traditional estimation of the variance of the 
bispectrum in the Gaussian limit [59H61I ]. i.e. Var[Sz 1 z 2 z 3 ] oc C l z 1 C'z 2 C’z 3 , with C; the angular power spectrum which 
is of the form Ci oc 1/1(1 + 1) in the range l < 20. On the other hand, if A 7 0. which physically means that the 
collapse of the wave function correlates different modes, then our expression (1561) . contains a function G(l) that is 
absent in the traditional approach. Furthermore, as we mentioned previously, in our approach there is in principle no 
restriction on l\ + I 2 + I 3 (which in the standard approach must be even a and represents an extra requirement for the 
nonvanishing of the bispectrum) as we are not assuming that the ensemble average of the bispectrum coincides with 
the observed one, this is, we do not basin our analysis in anything analogous to the standard claim that the 3-point 
correlation function of the temperature anisotropies must be rotational and translational invariant. This aspect of our 
analysis, represent a difference in the predictions that can be confronted with observations. In the next section, we 
will extend this and other discussions about the relation between the statistical aspects of the traditional bispectrum 
and the collapse bispectrum. 












IV. MAIN DIFFERENCES BETWEEN THE STANDARD AND THE COLLAPSE APPROACH 

REGARDING THE PRIMORDIAL BISPECTRUM 
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A. The standard approach to the CMB bispectrum 

We begin this section by giving a rather brief review of the conventional approach to the study of the primordial 
bispectrum, its amplitude and the usual arguments given to relate it with possible non-Gaussian features in the 
perturbations; extended reviews can be found in Refs. pH IdsUblj . 

According to the standard approach, if T (x, ij) is taken to be characterized by a Gaussian distribution all its 
statistical properties are codified in the two-point correlation func tion. Otherwise, one needs to consider higher 
order correlation functions, e.g. the three-point correlation function \k(:r, 77 )(5", rj). The Fourier transform is 

commonly referred as the bispectrum, defined by 


= (2tt) 3 (5(£i + k 2 + k 3 )By(ki,k 2 ,k 3 ). 


(58) 


The Dirac delta appearing in Eq. (1551) is said to indicate that the ensemble average 'S>{x,rj)'$!(y,r})'$!(z,ri) is invari¬ 
ant under spatial translations; in addition, the dependence only on the magnitudes k, appearing in the function 
By(ki,k 2 ,k 3 ), is tied to the rotational invariance of such ensemble. The Dirac delta constrains the three modes 
involved; this is, the modes must satisfy k± + k 2 + k 3 = 0, which is known as the triangle condition. Thus, according 
to the direct relation of the Newtonian potential with the source of the observed anisotropies, the assumption about 
the rotational and translational invariance of the ensemble and its impact on T(5?, ij) should be reflected 

in the invariance of the ensemble average of the temperature fluctuations ipi)^(0 2 , ^p 2 )^(d 3 , tp 3 ). However, 

we should once again emphasize that strictly speaking the discussion above and, thus, the average indicated by the 
overline in the latter expressions, refers to an average over an ensemble of universes and not averages over orientations 
in one universe. 

One of the first (and most popular ways) to parameterize the traditional non-Gaussianity phenomenologically is 
via the introduction of a non-linear correction to the linear Gaussian curvature perturbation [39l . |40| , 


= 'b g (x, V )+f^ 2 g (x, V ) - ^ 2 g (x,ri)\, (59) 


where 'k g (x , rj) denotes a linear Gaussian part of the perturbation and ^ g (x, rj) is the variance of the of the Gaussian 
part 03 The parameter fjjj is called the “local non-linear coupling parameter” and determines the “strength” of the 
primordial non-Gaussianity. This parametrization of non-Gaussianity is local in real space and therefore is called 
“local non-Gaussianity.” Using Eqs. (155)) and (1551) . the bispectrum of local non-Gaussianity may be derived: 


B*(fci, k 2 , k 3 ) — 2fjj'i J lPy(ki)Py(k 2 ) + Py{k 2 )Py(k 3 ) + Py(k 3 )Py(ki)\ 

_ 9 fioc f fcj k 2 k 3 \ 

NL {k\k 2 k 3 ) 2 \k 2 k 3 k\k 3 k\k 2 /’ 


(60) 


with P^(fc) the power spectrum of the Newtonian potential defined as SR (A:) (fc') = ( 27 t) 3 i 5(k + k')Py(k) and it is 

assumed to be of the form P^(k) = A^/c -3 , where Ay denotes the amplitude of the power spectrum. 

Expressing the coefficients ai m in terms of Tg as in Eq. (1551) and using Eq. (1551) . one can proceed to compute the 
CMB bispectrum, defined in Eq. (l40l) . this is 


16 Let us recall that the variance 'S/ 2 (x,rj) = dkk 2 P^(k,r]) diverges logarithmically for a power spectrum such that P\f,(k,rj) oc k~ 3 

unless one introduce an ad hoc cutoff for k. For a detailed discussion of this and other related issues see Ref. fl8ll and Appendix B of 
the same Ref. 
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(61) 


where in the last line, one performs the integral over the directions of the three vectors ki and use the exponential 
integral form for the delta function that appears in the bispectrum definition (1581) . The integral over the angular part 
of x, in the last line, is known as the Gaunt integral, i.e., 


/ d^Y limi (x)Y h m 2 (m 3 m 3 (x) 

I (2 li + l)(2l 2 + 1 )( 2 Z 3 + 1 ) / l\ l 2 h \ f h h h \ /pn', 

" V 4tt V 0 0 0 ) V m i m 2 m 3 ) ' 

The fact that the bispectrum B turns out to be proportional of the Gaunt integral, , implies that 

the values of l,m, corresponding to the non-vanishing components of the bispectrum, must satisfy “the triangle 
conditions” and also reflect the fact that in the definition of the bispectrum there is a Dirac delta [Eq. (1551) ] which 
guarantees the rotational invariance for the ensemble average 'p '■p. 

One can find a close form for the bispectrum by choosing to work in the Sach-Wolfe approximation, where the 
transfer function T(k) = 1; thus, substituting Eq. (1601) into Eq. (RJTl) and evaluating the remaining integrals one 
obtains (See Ref. |63j]): 
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\27tt 2 ) \h(h + i)h{h + 1) h{i2 +1 )h(k + 1) 


h{h + 1 ) ^3 (^3 + 1 ), 

Finally, using the definition of the angle-averaged-bispectrum [Eq. (l4ll) ]. one has 


(63) 


(2Zi + 1)(2^2 + 1)(2?3 + 1) f h h I 3 \ pioc ( 2A|, \ 

47t V° 0 0; JnL V27W 
1 1 1 
Will + ^)h{h + 1) h{h + l)l 3 {h + 1) h(li + l)l 3 (l 3 + 1) 

When working within the standard scenario for slow-roll inflation (and thus a “nearly” scale-invariant spectrum), 
Maldacena found that the estimate for the amplitude of non-Gaussianities of the local form, is ~ e [44j (this is 
assumed to be in the limit when k\ <C k 2 ~ k 3l i.e. in the so-called “squeezed” configuration). As seen from the 
previous discussion, in the conventional approach, the amplitude of the primordial bispectrum and the non-Gaussian 
statistics for the curvature perturbation are intrinsically related. 





B. Comparing the magnitude of the collapse and the traditional bispectrum 


The magnitude of \Bi 1 i 2 i 3 \ is obtained from Eq. flMt 
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(65) 


where we have used the estimate for the amplitude of the power spectrum for a single scalar field in the slow-roll 
scenario, given by Ay cs H 2 /(Mpt). 

On the other hand, the magnitude of the collapse biscpectrum is given by Eq. ©, which we will write again 
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( 66 ) 


Therefore we have two distinct theoretical predictions for the actual observed bispectrum \Bi x i 2 i 3 Actual obs [see Def. 
( 1571 ) ]. In the standard slow-roll inflationary scenario the prediction is given by |-B;p 2 z 3 | Eq. ( 1551 ) : meanwhile, by 
considering the collapse hypothesis the prediction is |Z?z 1 z 2 z 3 |m.l. Eq. ( 1551 ) . 

The first and most important difference is the fact that Eq. ( 1551 ) vanishes unless 7 ^ 0 , while no such “primordial 
non-Gaussianity” is required for a non-vanishing value of Eq. ( 1551 ) . 

The second difference is that the shape of the bispectrum, i.e. its dependence on l, is not the same; |m.l. 

scales roughly as ~ [(1 + A G(l))/l(l + l )] 3 / 2 while l-Bzp^l as ~ (21 + l) 3 / 2 /[l 3 (l + l) 2 ]. Unfortunately, the existing 
analysis of the observational data do not focus the exact shape of the biscpectrum, but rather on a generic measure 
of its amplitude. However, the reported observational amplitude of the bispectrum, which in the standard picture of 
slow-roll inflation corresponds to the non-linear parameter depends on the expected shape that emerges from 
the theoretical estimates of bispectrum [49|, l5(ij ]. In other words, in the standard approach, in order to obtain an 
estimate for the amplitude of the bispectrum from the observational data, one requires a theoretical motivated shape 
for the bispectrum. The observational data and theory are strongly interdependent. In this way, and relying on 
such theoretical considerations, the latest results from Planck mission [62| lead to an estimate for the amplitudes of 
the bispectrum for the local, equilateral, and orthogonal models given by /^l = 2.7 ± 5.8, /^ ral = —42 ± 75, and 
/NL ho = ~25 ± 39 ( 68 % CL statistical). On the other hand, none of the previously considered shapes, namely the 
local, equilateral or orthogonal, correspond to the one given by the collapse bispectrum. 


C. Comparison with the variance of the traditional bispectrum 


Furthermore, one can compute the cosmic variance of the traditional bispectrum, 
Gaussianity. This is, in the Gaussian limit, the variance of the standard bispectrum is [51 


in the case of weak non- 
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Meanwhile rewriting Eq. © yields 
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( 68 ) 


Comparing Eqs. © and ©, we note some similarities but also fundamental differences: First, if h ^ I 2 ^ I 3 
and the sum h + I 2 + h = odd, our theoretical prediction is |£>ZiZ 2 z 3 Im l = D^D^D^, meanwhile in the standard 
approach the theoretical prediction would be \Bi 1 i 2 i 3 \ 2 = 0. This is because, the traditional approach makes the 
additional assumption that the 3-point correlation function must be invariant under spatial translations and rotations; 
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in turn, in our approach, since we do not identify theoretical averages with observational quantities, the selection rule 
h + h + h = even, does not apply. Moreover, if l\ + l 2 + I 3 = e ven and A = 0, then the variance of the traditional 
bispectrum and our prediction matches exactly, this is, |S ; 1 ; 2 ; 3 | 2 = |B;i; 2 ; 3 |ml ^ ChCi 2 C; 3 . Note, however, that in 
general A 7 ^ 0 since this parameter characterize the corre lation between modes that is induced by the collapse of the 
inflaton’s wave function. Therefore, generically |-B / 1 ; 2 ; 3 | 2 7 ^ I^iwJml i moreover, there is the issue of the parity of 
h + h + h, mentioned above , all of which lead to the difference in the estimations of the two approaches reflected in 
the two expressions Eqs. (EZl) and (1551) . 

We must of course acknowledge that adopting a pragmatic point of view would lead one to argue that within 
the standard approach, and in the absence of non-gaussianities, the expected bispectrum would be zero up to the 
cosmic variance that can be seen in the observational data. In such case, our prediction, might be reinterpreted as 
an “improved estimate” of the traditional cosmic variance associated to the bispectrum, and in fact, the effects for 
introducing the collapse of the wave function, encoded in the function G(l), can be sought within the usual bispectrum 
cosmic variance. 

However, from the conceptual point of view, there are several differences between ours and the standard approach. 
Within the collapse hypothesis perspective, the observed \Bi x i 2 i 3 | Actual obs is clearly identified as corresponding to 
just one particular realization of a random quantum process (the self-induced collapse of the wave function). Since 
we do not have access to other realizations (i.e. we do not have observational access to other universes) we cannot 
say anything conclusive as to whether the underlying PDF is Gaussian or not. In other words, by measuring a 
non-vanishing |Z3; 1 ; 2 ; 3 | Actual obs in our own universe, which corresponds to a single realization of the physical process, 
cannot be taken as indication that the ensemble average , is also non-vanishing and would consequently 

proving non-Gaussianity statistics for the ensemble. In fact, just as we do not expect the actual value of the one 
available realization of a; m (for a fixed value of l and to) to vanish identically, even if somehow the ensemble average 
of such quantity (something to which we would have no access even if an ensemble of universes did exist) vanishes, we 
should not expect the single realization of a bispectrum, corresponding to our observations of the universe, to vanish 
identically, even if the average value over an ensemble of universes would vanish. 


V. CONCLUSIONS 

We have presented a detailed discussion on the manner, in which the study of essential statistical features on the 
CMB spectrum, must be studied in the context of the collapse models. We have shown the important differences that 
arise between the analysis in this approach and those tied to the standard approaches. We have seen, for instance, that 
the collapse models lead to explicit expressions for the quantities that are rather directly observable, the ai m [see Eq. 
m, which have no counterpart in the usual analyses. Among other advantages, the expression for the coefficients 
ai m [Eq- (1271) ] exhibits directly the source of the randomness involved, aspects that in the standard approach can 
only be discussed heuristically (simply because there, the random variables are not clearly identified and named as 
in this approach^). We have seen that this leads to some differences in the analysis of the higher order statistical 
features, such as the angular power spectrum [Eq. (1391) ] and the bispectrum [Eq. (1 55 f) ]. Within the standard approach, 
the bispectrum is usually associated with non-Gaussianities, meanwhile in our approach, it is a characterization of 
the limited observational data associated with the single realization of an explicitly described random process. On 
the other hand, the usage of the bispectrum to infer statistical aspects of an hypothetical ensemble of universes, 
characterized by a multiplicity of realizations of primordial perturbations, can be seen as a involving rather delicate 
extrapolations. 

We have seen that in the collapse scheme, one is lead to expect generically a non-vanishing bispectrum, even if 
there are no primordial non-Gaussianities (i.e. there are no assumptions regarding a non-vanishing quantum-tree 
point function, and in the statistical analysis of Sec. Mil all statistical correlations are encoded in the two variable 
correlation functions of the fiducial ensemble of universes, which are taken as standard). On the other hand, within 
the traditional approach, if one assumes that the primordial curvature perturbation is Gaussian, then the bispectrum 
must be zero up to its cosmic variance Eq. E3- Furthermore, after comparing our prediction for the observed 
bispectrum and the traditional variance associated to the bispectrum, we found that in general these two expressions 
are not equal (although for certain particular cases they do coincide). In fact, some of the differences arise directly 
from the assumptions regarding the self-induced collapse of the wave function: First, there is the issue of the parity 
of the bispectrum, in the traditional approach l± + l 2 + I 3 must be even, for a non-vanishing value of the bispectrum 
in order to be in agreement with a statistical homogeneity and isotropy of the ensemble of temperature anisotropies. 
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Researchers involved in direct comparison of actual empirical data with theoretical estimates apparently make use of similar random 
parameters in performing the simulations used in the process of confrontation of theory and observations f67t| 
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In our model, there is no constraint on the parity of l± -f 12 + 13 as we do not associate in any direct way the average 
over an ensemble of universes with the actual observed value of the bispectrum, which corresponds to just a particular 
element of the hypothetical ensemble. Second, the collapse of the wave function induces small correlations between 
different modes characterizing the inflaton (where the “smallness,” is parameterized by A), which in turn, imprints a 
particular signature in both, the angular spectrum (see Fig. 2) and the bispectrum that is encoded in the function 
G(l) [see Eqs. (l34|) and (j68|) ]. Therefore, even if one takes a pragmatical point of view, and argue that what we 
have actually computed is “just” the variance of the bispectrum, for the case when the primordial perturbations are 
Gaussian, our result differs in various mathematical aspects from that obtained in the standard approach. These 
differences (mainly those associated with the parity and the additional function G(l)) indicate that our prediction 
must at least be viewed as a novel one and we think it constitutes an improved estimate of the bispectrum cosmic 
variance. This makes the analysis highly predictive, and by the same token, highly susceptible to falsification by 
observational data, i.e. by using the data analysis, say from Planck collaboration, one could search for a signal that 
matches our prediction; in particular, the function G(l) mostly affects the lowest multipole angular scales l < 20 and 
then it decreases its effect on large multipoles in a nearly exponential fashion (see Figs. 1,2). 

Nevertheless, before even proceeding to do this, one would need to reevaluate the predicted angular spectrum 
and bispectrum by taking into account the effect of the transfer functions that we have ignored in performing the 
calculation leading to the expressions (1391) and (l56l) . That is, in such calculation, we replaced the transfer functions by 
the number 1 so that the integrals could be evaluated in closed form. Thus, the actual comparison of the prediction 
of this model with data will require the reintroduction of the transfer functions in the evaluation. Carrying this out 
would involve numerical calculations in order to perform the desired integrals that will result in the specific form 
of the Bi x i 2 i 3 , which will be suitable for comparison with observations. We plan to carry this analysis in the near 
future and to obtain the data to contrast with the model’s predictions. We would view a reasonable match between 
observations and the model as a strong indication we are on the right track. 

We end up by noting that there are a couple of aspects of our analysis that have suggestive counterparts in 
observations that have from time to time caught the attention of some researchers in the field: 

1. There is a claim of an abnormal deficit in power in the spectrum at low values of Vs [68l [69| which we think can 
be naturally explained within our approach by a negative value of A in equation (15^1) . 

2. There has been some attention given to an anomalous high degree of alignment occurring in the orientation of 
the low l multipoles (sometimes called the “axis of evil” [70]) and it seems again that a high correlation between 
the modes k and 2k that we have explained is a natural result of the collapse mechanism as viewed within the 
SSC formal analysis, would have the basic features required to account for that. 

These and other points would of course deserve further inquiry that falls outside the scope of the present manuscript. 


ACKNOWLEDGMENTS 

We wish to thank Prof. Matt Visser and Prof. Arthur Kosowsky for useful suggestions that improved the presenta¬ 
tion of the results obtained in this paper. The work of GL is supported by postdoctoral grant from Consejo Nacional de 
Investigaciones Cientificas y Tecnicas, Argentina. The work of DS was supported in part by the CONACYT-Mexico 
Grant No. 101712 and UNAM-PAPIIT Grant No. IN107412. 


[1] V. F. Mukhanov and G. V. Chibisov, JETP Lett. 33, 532 (1981). 

[2] S. W. Hawking, Phys. Lett. B 115, 295 (1982). 

[3] A. A. Starobinsky, Phys. Lett. B 117, 175 (1982). 

[4] A. H. Guth and S.-Y. Pi, Phys. Rev. Lett. 49, 1110 (1982). 

[5] J. J. Halliwell and S. W. Hawking, Phys. Rev. D 31, 1777 (1985). 

[6] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 (1992). 

[7] G. Hinshaw et al. [WMAP Collaboration], Astrophys. J. Suppl. 208, 19 (2013). 

[8] C. P. Ahn et al. [SDSS Collaboration], Astrophys. J. Suppl. 203, 21 (2012). 

[9] P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 571, (2014) A22. 

[10] A. Albrecht, P. Ferreira, M. Joyce and T. Prokopec, Phys. Rev. D 50 (1994) 4807 [astro-ph/9303001]. 

[11] D. Polarski And A. A. Starobinsky, Class. Quant. Grav. 13 (1996) 377 [gr-qc/9504030]. 

[12] D. Sudarsky, Int. J. Mod. Phys. D 20 (2011) 509 [arXiv:0906.0315]. 

[13] A. Perez, H. Sahlmann and D. Sudarsky, Class. Quantum Grav. 23, 2317 (2006). 



24 


[14] T. Padmanabhan, Structure formation in the universe (Cambridge University Press, New York, 1993). 

[15] V. F. Muklranov, Physical foundations of cosmology (Cambridge University Press, New York, 2005). 

[16] S. Weinberg, Cosmology (Oxford University Press, New York, 2008). 

[17] D. H. Lyth and A. R. Liddle, The primordial density perturbation: cosmology, inflation and the origin of structure (Cam¬ 
bridge University Press, New York, 2009). 

[18] G. Leon, S. J. Landau and D. Sudarsky, Phys. Rev. D 88, 023526 (2013). 

[19] P. M. Pearle, arXiv:0710.0567. 

[20] D. Sudarsky, J. Phys. Conf. Ser. 68, 012029 (2007). 

[21] D. Sudarsky, J. Phys. Conf. Ser. 67, 012054 (2007). 

[22] D. Sudarsky, Int. J. Mod. Phys. D 20, 821 (2011). 

[23] L. Diosi, Phys. Lett. A 105, 199-202, (1984). 

[24] P. M. Pearle, Phys. Rev. D 13, 857-868, (1976). 

[25] P. M. Pearle, Phys. Rev. A 39, 2277 (1989); G. C. Ghirardi, P. M. Pearle and A. Rimini, Phys. Rev. A 42, 78 (1990). 

[26] G. C. Ghirardi, A. Rimini and T. Weber, Phys. Rev. D 34, 470 (1986). 

[27] R. Penrose R, The Emperor’s New Mind, (Oxford University Press, New York, 1989), p. 480; R. Penrose, Gen. Relativ. 
Gravit. 28, 581 (1996). 

[28] L. Diosi, Phys. Lett. A 120, 377 (1987); L. Diosi, Phys. Lett. A 40, 1165 (1989). 

[29] S. Weinberg, Phys. Rev. A 85, 062116 (2012). 

[30] A. Bassi and G. C. Ghirardi, Phys. Rept. 379, 257 (2003). 

[31] P. Canate, P. Pearle and D. Sudarsky, arXiv:1211.3463. 

[32] J. Martin, V. Vennin and P. Peter, Phys. Rev. D 86, 103524 (2012). 

[33] S. Das, K. Lochan, S. Sahu and T. P. Singh, Phys. Rev. D 88, 085020 (2013). 

[34] A. De Unanue and D. Sudarsky, Pliys. Rev. D 78, 043510 (2008). 

[35] S. J. Landau, C. G. Scoccola and D. Sudarsky, Phys. Rev. D 85, 123001 (2012). 

[36] D. Larson, J. Dunkley, G. Hinshaw, E. Komatsu, M. R. Nolta, C. L. Bennett, B. Gold and M. Halpern et al., Astrophys. 
J. Suppl. 192, 16 (2011). 

[37] B. A. Reid, W. J. Percival, D. J. Eisenstein, L. Verde, D. N. Spergel, R. A. Skibba, N. A. Bahcall and T. Budavari et al., 
Mon. Not. Roy. Astron. Soc. 404, 60 (2010). 

[38] E. Komatsu, N. Afshordi, N. Bartolo, D. Baumann, J. R. Bond, E. I. Buchbinder, C. T. Byrnes and X. Chen et al, 
arXiv:0902.4759. 

[39] D. S. Salopek and J. R. Bond, Phys. Rev. D 42, 3936 (1990). 

[40] D. S. Salopek and J. R. Bond, Phys. Rev. D 43, 1005 (1991). 

[41] T. Falk, R. Rangarajan and M. Srednicki, Astrophys. J. 403, LI (1993). 

[42] A. Gangui, F. Lucchin, S. Matarrese and S. Mollerach, Astrophys. J. 430, 447 (1994). 

[43] V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, Nucl. Phys. B 667, 119 (2003). 

[44] J. M. Maldacena, JHEP 0305, 013 (2003). 

[45] N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402, 103 (2004). 

[46] A. Diez-Tejedor and D. Sudarsky, JCAP 1207, 045 (2012). 

[47] G. Leon and D. Sudarsky, SIGMA 8, 024 (2012). 

[48] E. Komatsu, Class. Quantum Grav. 27, 124010 (2010) 

[49] A. P. S. Yadav and B. D. Wandclt, Adv. Astron. 2010, 565248 (2010). 

[50] M. Liguori, E. Sefusatti, J.R. Fergusson and E.P.S. Shellard, Microwave Background and Large-Scale Structure, Adv. 
Astron. 2010, 980523 (2010). 

[51] E. Komatsu, Pli.D. thesis, Tohoku University, (2001). 

[52] G. Leon, A. De Unanue and D. Sudarsky, Class. Quantum Grav. 28, 155010 (2011). 

[53] A. Diez-Tejedor, G. Leon and D. Sudarsky, Gen. Relativ. Gravit. 44, 2965 (2012). 

[54] A. Valentini, Phys. Rev. D 82, 063513 (2010). 

[55] N. Pinto-Neto, G. Santos and W. Struyve, Phys. Rev. D 85, 083506 (2012). 

[56] C. J. Isham, In *Salamanca 1992, Proceedings, Integrable systems, quantum groups, and quantum field theories* 157-287, 
and London Imp. Coll. - ICTP-91-92-25 (92/08,rec.Nov.) 124 p [gr-qc/9210011]. 

[57] M. Fox, Quantum Optics An Introduction (Oxford Master Series in Physics, Oxford University Press, New York, 2006). 

[58] G. Leon and D. Sudarsky, Class. Quantum Grav. 27, 225017 (2010). 

[59] X. c. Luo, Astrophys. J. 427, L71 (1994) 

[60] A. Gangui and J. Martin, Mon. Not. Roy. Astron. Soc. 313, 323 (2000) 

[61] A. Gangui and J. Martin, Phys. Rev. D 62, 103004 (2000) 

[62] P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 571, A24 (2014) 

[63] J. R. Fergusson and E. P. S. Shellard, Pliys. Rev. D 80, 043510 (2009). 

[64] S. Weinberg, Phys. Rev. D 72, 043514 (2005). 

[65] X. Chen, Adv. Astron. 2010, 638979 (2010). 

[66] X. Chen, M. -x. Huang, S. Kachru and G. Shiu, JCAP 0701, 002 (2007). 

[67] S. Aiola, A. Kosowsky and B. Wang, Phys. Rev. D 91 (2015) 4, 043510. 

[68] C. J. Copi, D. Huterer, D. J. Schwarz and G. D. Starkman, Adv. Astron. 2010, (2010) 847541. 

[69] P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 571, (2014) A15. 

[70] K. Land and J. Magueijo, Phys. Rev. Lett. 95, (2005) 071301. 



